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Abstract. The form factor bootstrap approach allows to construct the space of local fields in the massive 
restricted sine-Gordon model. This space has to be isomorphic to that of the corresponding minimal model 
of conformal field theory. We describe the subspaces which correspond to the Verma modules of primary 
fields in terms of the commutative algebra of local integrals of motion and of a fermion (Neveu-Schwarz or 
Ramond depending on the particular primary field). The description of null- vectors relies on the relation 
between form factors and deformed hyper-elliptic integrals. The null-vectors correspond to the deformed 
exact forms and to the deformed Riemann bilinear identity. In the operator language, the null-vectors are 
created by the action of two operators Q (linear in the fermion) and C (quadratic in the fermion). We show 
that by factorizing out the null-vectors one gets the space of operators with the correct character. In the 
classical limit, using the operators Q and C we obtain a new, very compact, description of the KdV hierarchy. 
We also discuss a beautiful relation with the method of Whitham. 
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1 Introduction. 



In this article we present a synthesis of the ideas of the papers [|T| and j2|. In the first of these papers the 
space of fields for the sine-Gordon model (SG) was described in terms of the form factors previously obtained 
in the bootstrap approach ||. This description is based on rather special properties of form factors for the 
SG model. Namely, it uses the fact that the form factors were written in terms of deformed hyper-elliptic 
differentials, allowing deformations of all the nice properties of the usual hyper-elliptic differentials: the notion 
of deformed exact forms and of the deformed Riemann bilinear identity are available for them 0. Using 
these facts it has been shown in p| that the same number of local operators can be constructed in the generic 
case of the sine-Gordon model as at the free fermion point. The deformed exact forms and the deformed 
Riemann bilinear identity are necessary in order to reduce the space of fields to the proper size because in 
its original form factor description the space is too big. The description of (I) is basically independent of the 
coupling constant, but for rational coupling constant there is a possibility to find additional degenerations. 

The problem with the description of the space of fields obtained in |l| is due to the fact that it is difficult 
to compare it with the description coming from Conformal Field Theory (CFT) or from the classical theory. 
The latter two are closely connected because the Virasoro algebra can be considered as a quantization of the 
second Poisson structure of KdV. In the description of Q , it is even difficult to distinguish the descendents 
with respect to the two chiral Virasoro algebras. 

On the other hand in ]2| the semi-classical limit of the form factor formulae has been understood. This 
opens the possibility of identifying all the local operators by their classical analogues. Using this result we 
decided to try to construct the module of the descendents of the primary fields with respect to the chiral 
Virasoro algebra. The result of this study happened to be quite interesting. 

Let us formulate more precisely the problems discussed in this paper. The sine-Gordon model is described 
by the action: 

S = - J ({d^ip) 2 + m 2 {cos(2tp) - 1)) d 2 x 

where 7 is the coupling constant. In the quantum theory, the relevant coupling constant is: £ = -^zj- 

The sine-Gordon theory contains two subalgebras of local operators which, as operator algebras, are 
generated by exp{iip) and exp(— iip) respectively. We shall consider one of them, say the one generated by 
cxp(i(p). It is known that this subalgebra can be considered independently of the rest of the operators, as 
the operator algebra of the theory with the modified energy-momentum tensor: 

where a — tt^ This modification changes the trace of the energy-momentum tensor which is now: 
T™° d = m 2 exp(2i<^>). This modified energy-momentum tensor corresponds to the restricted sine-Gordon 
theory (RSG). For rational -, the RSG model describes the "Sji^j-perturbations of the minimal models of 
CFT. In this paper we consider only the RSG model. 

It is natural from the physical point of view of integrable perturbations || to expect that the space of 
fields for the perturbed model is the same as for its conformal limit. The latter consists of the primary fields 
and their descendents with respect to the two chiral Virasoro algebras. In this paper we shall consider the 
descendents with respect to one of these algebras, the possibility of considering the descendents with respect 
to the other one is explained in Subsection 2.2. The Verma module of the Virasoro algebra is generated by 
the action of the generators L_fc of the Virasoro algebra on the primary field. The irreducible representation 
corresponding to a given primary field is obtained by factorizing out the null- vectors |], |?| . 

On the other hand the very possibility of integrable deformations is due to the fact that there exists a 
commutative subalgebra of the universal enveloping algebra of the Virasoro algebra: the algebra of local 
integrals || ||, |). The local integrals /2fc-i have odd spins. Logically it must be possible to present the 
Verma module as a result of the action of the local integrals of motion on a smaller module: the quotient over 
the local integrals. The latter is isomorphic to the module created from the vacuum by the action of bosons 
of even spins J2/C ■ The important observation is that the form factor formulae allow us to give a description 
of this space. Certainly, such a realization of the Verma module is very useful for integrable applications, 
but there is a difficulty. It is not trivial to describe the null-vectors in this realization. This question, of the 
description of the null-vectors, is the main problem solved in this paper. 

Briefly, the result is as follows. It is useful to fermionize the bosons J2k by introducing Neveu-Schwarz 
(ip2k-i) or Ramond (ip2k) fermions, depending on which primary field we consider. As in [Q the null-vectors 
are due to the deformed exact forms and the deformed Riemann bilinear identity. In the fermionic language 
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these null- vectors are created by the action of two operators: a linear one in the fermion (Q) and a quadratic 
one (C). By calculating the characters we show that, after factorizing these null- vectors, we find exactly the 
same number of dcsccndcnts of a given primary field as in the corresponding irreducible representation of 
the Virasoro algebra. 

In Conformal Field Theory, the existence of null vectors provide differential equations for the correlation 
functions. The null vectors built out of Q and C in our approach will also lead to a system of equations for 
the correlation functions in the massive case. These equations will be similar to the hierarchies of equations 
in classical integrable systems. In spite of the fact that such an infinite set of equations seems a priori 
untractable, our hope is that, as in the classical case, large classes of solutions will be constructed. Hence 
we hope that these equations will be really useful for the study of correlation functions in the massive case. 

We also present a detailed analysis of the classical limit of our constructions. The description of the 
null- vectors in terms of Q and C implies in the classical limit, a new and very compact description of the 
KdV hierarchy. This new description is not the same as the description in terms of r-functions [fit)) , still 
the technics we use are close to those of the Kyoto school. Finally, we discuss an amazing analogy between 
the quantum theory and the results of the Whitham method for small perturbations around a given quasi- 
periodic classical solution of KdV. 

Let us make one general remark on the exposition. In this paper we restrict ourselves to the reflectionless 
case of the sine-Gordon model. This is done only in order to simplify the reading of the paper. All our 
conclusions are valid for a generic value of the coupling constant as well. For this reason, we chose to 
present the final formulae in the situation of a generic coupling constant. This leads, from time to time, to 
a contradictory situation. We hope that we shall be forgiven for this because if we wrote generic formulae 
everywhere the understanding of the paper would have been much more difficult. 

2 The space of fields. 

2.1 The description of the space of fields in the A, B- variables. 

At the reflectionless points (£ = ^, v = 1, 2, • • •) there is a wide class of local operators O for which the form 
factors in the RSG-model corresponding to a state with n-solitons and n-anti-solitons are given by 

fo(fa,fa, ■■■,02n) +■■■+ = 

n 2n 1 1 

=c»nc(A-/%)n n sinh ^.- A - 7rt ) ^-2^- i )-^ft) 

i<j i=lj =n +l K ' rl ' j 

X fo(Pl, 02, ■■■,&*)- .- + ...+ (1) 

The function without poles in the strip < Im (3 < 2n, satisfies £(— /3) = S(/3)£((3) and ((f3 — 2ixi) = 
£(— 0): the S-matrix S{(3) and the constant c are given in the Appendix A. The most essential part of the 
form factor is given by 



A)) L<g\Ai, • • • , A n \Bi, • • • , B 2n ) aT (2) 
where Bj = e l3j and 

i/f(A,B) = Y[(B- Aq- j ), with q = e™ lu 

As usual we define a — A 2l/ . Here and later if the range of integration is not specified the integral is taken 
around 0. Notice that the operator dependence of the form factors d3) only enters in fo- 

Different local operators O are defined by different functions Lq (Ai, • • • , A n \B\, ■ • ■ , i?2n)- These func- 
tions are symmetric polynomials of Ai,---,A n . For the primary operators $2fc — exp(2kiip) and their 
Virasoro descendants, Lq are symmetric Laurent polynomials of Bi, ■ ■ ■ , i?2n- For the primary operators 

^2fc+i = exp((2fc + l)i(p), they are symmetric Laurent polynomials of B\,- ■ ■ , i?2n multiplied by FJ BJ . 
Our definition of the fields $ m is related to the notations coming from CFT as follows: $ m corresponds 



fo(Pl, fa, • ■ • , fan) + ■■■+ = 

n 2n 



(2m) 



' J J i=l9=l i<j 
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to 5?ri,m+i]- The requirement of locality is guaranteed by the following simple recurrent relation for the 
polynomials Lq\ 

L@\Ai, A n \Bi, B 2n ) = 

w B 2n =-Bi, A n =±B x 

= -e ± L^ 1) (A 1) ■ • ■ , A n _i|S 2) • • ■ , B 2n _ x ) (3) 

where e = + or — respectively for the operators <i>2fc and their descendents, or for <&2fc+i and their descendents. 
In addition to the simple formula (Q) we have the requirement 

n In 

res An=00 ( Y[ Y[ 1>(Ai, Bj) J{{Aj - A 2 } ) L%\Ai, • • • , A n \B x , • • • , B 2n ) a - k ] = 0, k > n + 1 (4) 

yi=l J=l i<j 

which is true in particular if degA n {Lo) < 2v. We explain this conditions in Appendix A. We see that the 
restriction (|j) disappears only in the classical limit v — ► oo. It will be clear later that this class of local 
operators is not complete for the reason that the anzatz (^) is too restrictive. With this anzatz we obtain 
the complete set of operators only in the classical limit. However there is a possibility to define the form 
factors of local operators which correspond to polynomials satisfying the relation (^) without any restriction 
of the kind (Q). To do that for the reflectionless points one has to consider the coupling constant in generic 
position (in which case the formulae for the form factors are much more complicated [[|) and to perform 
carefully the limit £ = - + e, e — > 0. An example of such calculation for £ = 7r is given in 0. We shall 
give more comment about this procedure later, here we would like to emphasize that the local operator can 
be defined for any polynomial satisfying (|^) but its form factors are not necessarily given by the anzatz (^) . 
Physically the existence of local operators for the reflectionless case whose form factors are not given by the 
anzatz (|J) is related to existence of additional local conserved quantities which constitute the algebra sl(2). 
In spite of the fact that the form factors of the form (jl]) do not define all the operators they provide a good 
example for explaining the properties valid in generic case. 

The explicit form of the polynomials Lq for the primary operators is as follows 

n 2n 

42(^i. ■■■,A n \B 1 ,.-., B 2n ) = n at n b v 

i=i j=i 

In this paper we shall consider the Virasoro descendents of the primary fields. We shall restrict ourselves 
by considering only one chirality. Obviously, the locality relation (||) is not destroyed if we multiply the 
polynomial Lq{A\B) either by l2k-i(B) or by J2k{A\B) with 

hk-i{B) = Q^L! ) s 2k-i(B), fc=l,2,-.. (5) 

J 2 k(A\B) = s 2k {A) - ^s 2k (B), fc=l,2,.-. (6) 
Here and later we shall use the following definition: 
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The multiplication by I 2 k~i corresponds to the application of the local integrals of motion. The normalization 

-,, 2fc-l 

factor ^_/ 2 k-i is introduced for further convenience. Since the boost operator acts by dilatation on A and 
B, hk-i has spin (2k — 1) and J 2 k has spin 2k. 

The crucial assumption which we make is that the space of local fields descendents of the operator $ m 
is generated by the operators obtained from the generating function 



2/i 



C m (t,y\A\B) = exp(J2t2k-ihk-i(B)+y 2 kJ2k(A\B)j n^II B i ] ( 7 ) 
fc>i \i=i j=i J 

This is our main starting point. This assumption follows from the classical meaning of the variables A,5 |||. 
We shall give classical arguments later, and the connection with CFT is explained in the next subsection. 
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2.2 The relation with CFT. 



The RSG model coincides with the <J>[ 13 ]-perturbation of the minimal models of CFT. For the coupling 
constants which we consider the minimal models in question are not unitary, but this fact is of no importance 
for us. It is natural from a physical point of view, to conjecture that the number of local operators in the 
perturbed model is same as in its conformal limit. 

We first need to recall a few basic facts concerning the minimal conformal field theories. At the coupling 
constant £, the conformal limit of the RSG models have central charge 

c=l-6- ^ 



The conformal limit of the fields <fr m = e lmv is identified with the operators ^>[i jm +i] in the minimal models. 
They have conformal dimensions: 

_ m(m£ - 2tt) 
Am _ 4(£ + tt) 

The refiectionless points, £ = ^ with v = 1,2, correspond the non-unitary minimal models A^i^+i), 
which are degenerate cases. At these points one has to identify <I> m with <&2v-m- 

The Virasoro Verma module corresponding to the primary field $ m is generated by the vectors 

At £ = ^ their structure is the same as for generic value of £. This means that for m integer, the Verma 
module possesses only one submodule. The so-called basic null- vectors, which we shall denote by r m , are 
the generators of these submodules. They appear at level m + 1. The first few are given by: 



r = L_i$ , 

Ti = (L-a + KL?.!)*!, k=-1-" 



T 3 = (L_ 3 + Kli-li-2 + K 2 i 3 -i)$2, Kl = -2-^4,K2- (A+L) ~ 



tt + SC 2£(tt + 3£) 
etc ■ ■ ■ 

Other null-vectors, whose set form the submodule, are created from the basic ones as follows 

L-k x L-k 2 • • • L-k N r m 

As a consequence, the character of the irreducible Virasoro representation with highest weight <& m is: 

1 - p m+1 



Xm{p) = 



iwa-p*) 



There exists an alternative description of these modules which is more appropriate for our purposes. 
Indeed, as is well known the integrability of the $[ 13 ]-perturbation is related to the existence of a certain 
commutative subalgebra of the Virasoro universal enveloping algebra generated by elements hk-i, polyno- 
mial of degree k in L n , such that [L , hk-i] = (2fc — l)/2fe-i and 

[hk-uhi-i] =0, = 1,2, - - ■ 

The first few are given by: 

h = L-u 

h = 2 L_„_ 2 L„_i, 
etc 

This is the subalgebra of the local integrals of motion. The meaning of these operators and of the l2k-i(B) 
which we had above is the same, so we denote them by the same letters. 
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Of course acting only with this subalgebra on $ m does not generate the whole Verma module. But the 
left quotient of the Verma module by the ideal generated by the integrals of motion produces a space whose 
character is ]^[ - >0 (1 — p 23 ') -1 . This space can be thought of as generated by some J 2 k with [Lq, J 2 fc] — 2fcJ 2 fc. 

We expect that there is the following alternative way of generating the Virasoro module. Namely in a 
spirit similar to the Feigin-Fuchs construction |6| , we expect that there exists an appropriate completion of 
the subalgebra of the integrals of motion by elements J 2 fc such that (i) the Verma module is isomorphic, as 
a graded space, to the space generated by the vectors 

exp(^£ 2 fc-i^2fc-i) : exp(^y 2/c J 2/c ) : *m (8) 

k>l fc>l 

where the double dots refer to an appropriate normal ordering, and (ii) that for certain normalization of ijak— 1 
and certain choice of J 2 fc and their normal ordering the generating functions (ffl) and (||) provide different 
realizations of the same object. We shall call the description given by (0) the A, ^-representation of the 
Virasoro module. The existence of this alternative description of the Virasoro module is very important for 
integrable applications. 

The main problem for this /, J description arises from the non-trivial construction of the null-vectors. 
This problem has two aspects. First one has to construct the null-vectors T m in terms of I, J. Then one 
also has to construct the submodule associated to it. This is not as trivial as in the standard representation 
since, while acting on T m with the Jak— 1 still produces a null-vector, acting with the J 2 k does not necessarily 
leads to a null- vector. We shall show that the whole null- vector submodule can be described in the A, B- 
representation. The proof is based on rather delicate properties of the form-factor integrals, the main of 
them being the deformed Riemann bilinear identity Q. 

Let us discuss briefly the problem of the second chirality. In |2| we have explained that the formula ([!]) 
has to be understood as a result of light-cone quantization in which X- is considered as space and 2;+ as 
time. We must be able to consider the alternative possibility ( x + - space, x_ - time), and the results of 
quantization must coincide. Where exactly the choice of hamiltonian picture manifests itself in our formulae? 
Consider the formula (|[). As it is explained in the fact that is a polynomial in A4 corresponds to 

the choice of x~ as space direction while the multiplier Y[ a 7 Z corresponds semi-classically to the choice 

i=l 

of trajectories under x+-flow. These are the two ingredients which change when the hamiltonian picture is 
changed. Indeed, analyzing the results of we find the following alternative description of the form factors. 
The form factors are given by the formulae (Q) with f@ replaced by 

ho(Pi, 02, ■ ■ ■ , fan) 1 — h = 

„ „ n In n 2n 

= — — / dAx ■ •• / dA n 11 1[ tj}(Ai, Bj)Y[(A? A]) K%\A U - ■ • , A n \B x , ■ • • , B 2n ) J] a p +1 J] bf 

where bj = Bf, K { £ ] is a polynomial in A t 1 and the replacement of Yi a i 1 by Yi a i 1+1 Y\^j 2 corresponds 
to the change of ^^-trajectories to a;_-trajectories. In particular for the primary fields we have 



In 

K%1 (Ax, ■ ■ ■ , A n \B u • • • , B 2n ) - J] A- m [] B- 

i=i i=i 



The descendents are obtained by multiplying by I-(2k-i)(B) and J-2k(A\B). The consistency of 
the two pictures requires that for the primary fields they give the same result: 

f* m (Pi,fa #*,)-...-+...+ - h^ m (Px,p 2 , • ■ • , 3-m) - 1 1 (9) 

This is a complicated identity which nevertheless can be proven. We do not present the proof here because it 
goes beyond the scope of this paper. It should be said, however, that the proof is based on the same technics 
as used below (deformed Riemann bilinear identity etc) . Using the equivalence of the two representations of 
the form factors we can consider the descendents with respect to and also mixed L^k descendents. 
A formula similar to (Q) holds for any coupling constant. However, at the reflectionless points there is another 
consequence of (Bl): it also shows that the operators <I> m and 3>2i/-m are identified as it should be. 
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3 Null- vectors. 



3.1 The null-polynomials. 

Null-vectors correspond to operators with all the form factors vanishing. Consider the integral 

r, r, n In n 

—— dA l --- dA n ]l]l^(A i ,B j )]l(A 2 i -A])L^\A u ---,A n \B 1 ,...,B 2n )l[a^ (10) 

V / J i—l j—l i<j i—l 

Instead of Lq \ we shall often use the anti-symmetric polynomials Mq': 

A4 n) (Ai, • ■ • , A n \B u - ■ • , B 2n ) = YliAj - A]) L%\A 1: ■ ■ ■ , A„| J B 1 , • • • , B 2n ) 

i<j 

(n) 

The dependence on B±, ■ ■ ■ , B 2n in the polynomials Mq will often be omitted. 

There are several reasons why this integral can vanish. Some of them depend on a particular value of the 
coupling constant or on a particular number of solitons. We should not consider these occasional situations. 
There are three general reasons for the vanishing of the integral, let us present them. 

1. Residue. The integral ( |To| ) vanishes if vanishes the residue with respect to A n at the point A n = oo 
of the expression 

In 

H^(A, u B 1 )a- n M^ ) (A 1 ,---,A n ) 

3=1 

Of course the distinction of the variable A n is of no importance because Mq (Ax, ■ • ■ , A n ) is anti-symmetric. 

2. "Exact forms." The integral ( |To| ) vanishes if Mg' (A%, ■ ■ ■ , A n ) happens to be an "exact form". 
Namely, if it can be written as: 

m£\a u ■ ■ ■ , A n ) = ^(-l) fc M(^, • • • ,Tk, ■ ■ ■ , A.) (Q(A k )P(A k ) - qQ(qA k )P(-A k )) , (11) 

k 

with 

P(A) = H(B j+ A) 

3=1 

for some anti-symmetric polynomial M(A\, • • • , j4„_i). Here and later A k means that A k is omitted. This is 
a direct consequence of the functional equation satisfied by ip(A, B) (See eq. (pgj) in Appendix B). For Q(A) 
one can take in principle any Laurent polynomial, but since we want M^ 1 ' to be a polynomial the degree of 
Q(A) has to be greater or equal — 1. 

3. Deformed Riemann bilinear relation. The integral ( |To| ) vanishes if 

AfW(Ai, • ■ • , A n ) = 53(-l)<+^M(Ai, • • • ,A h ■ ■ ■ ,Aj, ■ ■ ■ A n )C(A t ,A 3 ) 
where M(Ai, • • • , A n ^ 2 ) is an anti-symmetric polynomial of n — 2 variables, and C(A\, A 2 ) is given by 

C(Al > M) = T[A 2 {a^T^ {P{Ai)P{A2) P (- A ^ p (- A *)) + [P{-A l )P{A 2 ) P(A 1 )P(-A 2 ))| (12) 

This property needs some comments. For the case of generic coupling constant its proof is rather complicated. 
It is a consequence of the so called deformed Riemann bilinear identity Q . The name is due to the fact 
that in the limit £ — > oo the deformed Riemann bilinear identity happens to be the same as the Riemann 
bilinear identity for hyper- elliptic integrals. The formula for C(Ai,A 2 ) given in differs from ( |l2|) by 
simple "exact forms". Notice that the formula for C(Ai,A 2 ) does not depend on the coupling constant. For 
the reflect ionless case a very simple proof is available which is given in Appendix B. 

In order to apply these restrictions to the description of the null-vectors we have to make some prepa- 
rations. The expression for C(A\, A%) given above is economic in a sense that, as a polynomial of A\,A 2 it 
has degree 2n — 1, but it is not appropriate for our goals because it mixes odd and even degrees of A\,A 2 
while the descendents of primary fields contain only odd or only even polynomials. So, by adding an "exact 
forms" we want to replace C(A±, A 2 ) by equivalent expressions of higher degrees which contain only odd or 
only even degrees. 
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Proposition 1. The following equivalent forms of C(Ai, A 2 ) exist: 

C(A 1 ,A 2 ) ~ C e (A 1 ,A 2 ) ~ C (Ai, A 2 ), (13) 

here and Jater ~ means equivalence up to "exact forms". The formulae for C e (A%, A2) and C (Ai, A2) are 
as follows: 

° MU M) = J d ° 2dDl Te UJ \(Dl-Al) { Di-Al) - { Dl-Al) { Dl-Al)) ^ 

|D 2 |>|£>i| 
|D 1 |>|A 1 |,| J 4 2 | 



where 

1 9 „2fc-l 1 + 9 



1 \ 1 — 9 " 2fe-l V""* 1 

Te W - 2^ 1 1 2fc-i x ~z^7~_ 2fc 

fc=l H k=l H 



x 2k 



and 



C (A 1 ,A 2 ) = A 1 A 2 J dD 2 dD\T 



DA f P(D 1 )P(D 2 ) P(D 1 )P(D 2 ) 



|D 2 l>|D 1 | 
Dll>|A 1 |,|A 2 | 



D 2 J \(D(-A{)(Di~Ai) (D 2 - A{)(D{ - A 2 ) 



where 

" l-9 2fc 2fe ^l + 9 2fc - 1 



The proof of this proposition is given in Appendix B. 

The functions T e (x) and t {x) are not well defined when q r — 1 because certain denominators vanish. 
However the formula ([l3]) in which the LHS is independent of q implies that in the dangerous places we 
always find "exact forms" . So, for our applications these singularities are harmless. We shall comment more 
on this point later. 

3.2 The fermionization. 

The descendents of the local operators are created by I 2 k-i and J 2 k- This generates a bosonic Fock space. It 
is very convenient to fermionize J 2k . Let us introduce Neveu-Schwarz and Ramond fermions: %j) 2 k-\, V'lfe— 1 
and tp 2 k , tp2k ■ The commutation relations are as follows 

We prefer to follow the notations from jl(| than those coming from CFT, the reader used to CFT language 
has to replace ipm by ip- m - 

The vacuum vectors for the spaces with different charges are defined as follows. In the Neveu-Schwarz 
sector we have: 

V>2fc-i |2to — 1} = 0, for k > m, ip2k-i\^ m — 1) = 0, for k < m; 
(2m — l\ip2k—l = 0, for k < m, (2m — l|"02fc-i = 0, for k > to 
For the Ramond sector we have: 

^2fe|2m) = 0, for k > to, 1 2m) = 0, for k < to; 

(2TO|-02fc = 0, for k < to (2m\ip 2k — 0, for k > to 

We shall never mix the Neveu-Schwarz and Ramond sectors. The spaces spanned by the right action of an 
equal number of ^'s and V>*'s on the vector (p| will be called H*. The right action of ip sends H* to H* +2 . 
It is useful to think of the vector (p\ as a semi-infinite product 

(p| = • • • 1pp-i^p-2lpp 

Let us introduce generating functions for the fermions. The operators ip(A),ip* (A) are defined for the 
Neveu-Schwarz and the Ramond sectors respectively as follows 

00 00 

tl>(A) = £ A- 2fc +V 2fc -i, W)= E ^V 2 Vi; 

k—~oo k— — oo 

00 00 

k— — oo k— — oo 
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We shall use the decomposition of ip(A), ip*(A) into the regular and singular parts (at zero): 

i/>(A) = i>(A) reg + i>(A) sing , ip*(A) = f (A) re9 + ^*{A) smg 

where tp(A) reg and ip*(A) reg contain all the terms with non-negative degrees of A. 
Let us introduce the bosonic commuting operators h-2k for k > 1: 

oo 

h-2k = /J V'2j-iV'2fe+2j-i f or Neveu — Schwarz sector, 

j=— oo 
oo 

/i_2fc = J]] V , 2jV , 2fc+2j / or Ramond sector, 

i=— oo 

They satisfy the commutation relations: [h^2k, h*_ 2l ] = —kSk,i- We also have the following commutation 
relations between the fermions and the bosons: 

^)fc_ 2fe = (/i_ 2fc - A" 2fc )^(A) (15) 

The bosonic generating function C m (t, y\A\B) for the descendents of the operators can be rewritten 



as: 



(t, y\A\B) = exp \^t 2 k-ihk-i{B)\ (m-l|exp(^y 2fe /i* 2fe ) C m {A\B)\m-l) (16) 



where 



,fc>i / fc>i 



C m (A\B) = expl-J2h- 2k J 2k (A,B)) 1^11 B J 
\ fc>i / * J 



In other words, to have a particular descendent one has to take in the expression 

exp Ij^^k-ihk-iiB) C m {A\B)\m - 1) 



, fc>i 



the coefficient in front of some monomial in i 2 fc-i and to calculate the matrix element with some vector from 
the fermionic space H* n _ 1 . 

We can replace this bosonic expression by a fermionic one. Recall that as a direct result of the boson- 
fermion correspondence one has: 



i<j \ k>l i=l J 



- 1) = r(Ai) ■ ■ ■ r(A n )\m 2n A] 



m+2n — 1 



where the fermions are Neveu-Schwarz or Ramond ones depending on the parity of to. Using this fact the 
formula ( |l6|) can be rewritten for any n as follows 



Y[(Aj-A 2 3 ) £ m (A\B)\m-l}= g(B) ip*(Ai)--- t(i*(A n )\m -2n-l) \\ A' \[ H 

i<j i j 

where 

^ ' -2 k s 2k (B) I (III 



Vfc>i 



The occurrence of this operator is similar to that of the spin field in the description |T^] of CFT on hyper- 
elliptic curves. 

Let us now concentrate on the operator $o = 1 and its descendents. It means that we are working with 
the H^_ 1 subspace of the Neveu-Schwarz sector. The polynomials corresponding to the operators in this 
sector are even in A;. Let us describe the null- vectors due to the restrictions 1, 2, 3, from the previous 
subsection. To do that we shall need the results of the following three propositions. 
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Proposition 2. The set of polynomials of the form 

Mg\Au — ,A n ) = ^(-1)^'M(A 1; ■ ...A, .....A,.... A n )C e (Ai, Aj) (18) 

i<j 

coincides up to "exact forms" with the set of the matrix elements: 

<*_ 5 |C r(A 1 )---r(A n )\-2n-l)l[A 2 i n -\ V tf_ 5 e H*_ 5 (19) 

where 

C= J dD 2 J dD l P{D l )P{D 2 )Df n - 1 D^ 2n - 1 r e (J^j ^{D{)^{D 2 ) 

\D 2 \>\D t \ 

Proof. 

The vectors from the space H* can be written as 

where k p > ■ ■ ■ > k\ > N + 1. We shall call N the depth of (vE^. There are three possibilities for the matrix 
element (|l9| ) to differ from zero: 

1. The depth of is greater than — 2n — 1 

2. The vector (\P-5I is obtained from a vector whose depth is greater than — 2n — 1 by application 
of ip-2p-i'4 > -2q~i w hh q > p > n (i.e. there are two holes below — 2n — 1). 

3. The vector (^-sl is obtained from a vector (^-sl whose depth is greater than — 2n — 1 by application 
of if)*_2p-i with p> n (i.e. there is one hole below — 2n — 1). 

In the first case using the formula 

(-2n-l\ip(D)r(A)\-2n-l)=(^) n f °\^ , \D\ > \A\ 



A J {D 2 - A 2 ) ' 



and (14) one find 

(*_ 5 |C piAt) ■ ■ ■ r(A n )\ -2n- 1} Y[A 2n l = 

= ....A, .....A,.... A n )C e {A,A 3 ) 

i<j 

where 

M(A U ■ ■ • , A n _ 2 ) = <*_ 5 |^*(^i) • • • r{A n -2)\ -In- 1) ]jA 2n 1 

In the second case it is necessary that in the expression -llV^p-iV^g-i ^ the two holes below — 2n — 1 
are annihilated by C, the result is 

0M J dD 2 JdD 1 P(D 1 )P(D 2 )D- 2n+2p D~ 2n+2g T e (^j=0 

\D 2 \>\Di\ 

because the integrand is a regular function of D\ for p > n. 

In the third case it is necessary that in the expression {^-3 \ip*_ 2p -\ C the hole below — 2n— 1 is annihilated 
by C, the result is 

<*- 8 | J dD 2 J dD 1 P(D 1 )P(D2)D^ 2n - 1 D- 2n+2p r e (^-)ij(D 1 



D 



|£) 2 |>|Di| 



where the pairing of ip(D\) and ip*_ 2p _i is not considered because it produces zero for the same reason as 
above. In the matrix element we shall have the polynomials 



(-2n - 1| J dD 2 J dD 1 P(D 1 )P(D 2 )D^ 2n - 1 D- 2n+2p T e (^j 



r/>(D 1 W{A j )\-2n-l)A 



2n-l 



\D 2 \>\D 1 \ 



J dD 2 J dD 1 P{D 1 )P(D2)^^D- 2n+2p T e {^j = R 2n+2p (A,) 



|D 2 |>|-Di| 
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The polynomial R 2n +2p(A) is an even polynomial of degree 2n + 2p. Let us show that it is an "exact form" . 
From the calculations of Appendix B we have 



R 2 „+2 P (A) = J dD 2 J dD 1 P{D 1 )P(D 2 ) D 2 l _ A 2 D 

\D-z\>\Di\ ' 

^ A J dD 2 J dD 1 P(D 1 )P(D 2 ) D ^ 1 _ A2 a 



Ei 

D 2 



2n+2p 

2 



D 1 + D 



|0 2 |>|-Di| 



A J dD 2 J "^(W)^^ (^) 



(20) 



|Di|>|-D2 



where we have changed the integral over \D 2 \ > \Di\ by the integral over \Di\ > \D 2 \ because the residue at 
D 2 = —D\ gives the integral over D\ of an even function, such integral equals zero. The last integral in ( pc| ) 
vanishes because the integrand is a regular function of D 2 for p > n. Let us emphasize that our construction 
is self-consistent because for every n the polynomials of too high degree (greater than An — 2) are "exact 
forms" . 

Thus we have non-trivial matrix elements only in the first case which obviously exhausts the polynomial 
of the kind (£§). ■ 

Consider now the restriction 2 of the previous subsection. It is easy to figure out that there is only one 
uniform way to write for all n polynomials of the type (0) which are even in all variables Aj. Namely: 



Mg\A u ■ ■ ■ , A n ) = £(-l) fe MK, ■ • • , Al ■ ■ ■ , Al) (P(A k ) P(-A k )) Af (21) 
k 

where M(A\ , is arbitrary anti-symmetric polynomial. 

The following two simple propositions are given without proof. 

Proposition 3. Tiie set of polynomials ( pi} ) coincides with the set of matrix elements: 
(*_ 3 |Q r{Ax) ■ ■ ■ r{A n )\ - 2n - 1) I] A 2 "" 1 V*_ 3 e H* 3 

where 

Q = J dDD- 2n - 1 P(D)ijj(D) 

Proposition 4. The set of polynomials M^(Ai, • • • , A n ) such that 

res An=oc J] i/>(A n , B )a^ n M^\A ll - ■ ■ , A n ) = 

3 

(we hope that the same letter ip used for the function ip(A, Bj) and for the fermion is not confusing) coincides 
with the set of matrix elements: 

(* 1 \Q ] r{A 1 )---r{A n )\ -2n-l)l[A 2 l n - 1 V&a E H* 

where 

& =res A=00 l[^(A,B J )a- n J ^ (D)^E-^D 2n dD 

\D\>\A\ 

Let us apply these results to the description of null- vectors. We need to introduce the following notations: 

jik-i- 



k>l k>l v H 7 



2k 

k>l 
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Obviously 

P(D) = D ^ e X(D)-Y{D) (24) 

The null- vectors will be produced by acting with some operators C, Q and Q) on jCq(A\B) \ — 1). In view of 
the bosonization formulae, these operators are obtained from C, Q and QJ by conjugation with g(B): 

C g(B) = g(B) C, Q g(B) = g(B) Q, Qt g(B) = g(B) & 

The formulae for C and Q are given in the following two propositions: 

Proposition 2'. From the Proposition 2 we find that the null-vectors due to the deformed Riemann bilinear 
identity are of the form: 



exp(J2t2k-ihk-i(B)) (*_ 5 | C exp{-^2^h- 2k J 2k {A,B))\ - 1) 

k>l k>l 

where 

C= J / ^e^e^Ve (g) WMlh) (25) 

\D 2 \>\D 1 \ 

Proposition 3'. From the Proposition 3 one gets the following set of null-vectors due to the "exact forms": 
exp{J2t2k-il2k-i{B)) (9- 3 \Qexp{-^2-h-2kJ2k(AS))\-l) 



k 

fc>l k>l 



where 



Q = J d -§e^ D H{D) (26) 

Notice a very important feature in these formulae: The operators Q and C are independent of n. 
These propositions are direct consequences of the previous ones and of the following conjugation property 
of the fermions: 

1>(D) g(B) = g(B) ip(D) e~ Y ^\ ip{p) g(B) = g(B) e Y ^ 

Before dealing with let us discuss the operator C, Q in more details. It will be convenient to rewrite 
them in terms of another set of fermions ip and tp' . To understand the purpose of introducing a new basis 
for the fermions consider the formula (p5|). In this formula the fermion tp(D2) can be replaced by its regular 
part Tp{D2)reg because other multipliers in the integrand contain only negative powers of Z?2- That is why 
C can be rewritten in the form 



C = I ^-Tp(D) smg1 p(D) re g (27) 



where the modified fermion -0 is defined as follows 

4?{D) reg = xp(D) reg , 4>{D) slng = Ui> (D), 
with U the following operator 



Uf (D) = 



^£l p X{D 1 )X{D) T f£l 



D 1 \D 

\D\>\Di\ 



f(Di) 



odd 



where [• • -] dd means that only odd degrees of the expression with respect to D are taken because only those 
contribute to the integral (^7|). It is quite obvious that this transformation is triangular, namely 

/l_ q 2k-l\ 

tp2k-i = i i + q 2k-i J V^fc-i + (terms with ip 2 i-i,l < k), k>l 
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Altogether we can write ip(D) = Uip (D) where U is triangular. 

Introduce the fermions ip' satisfying canonical commutation relations with ip: 



tf(D)=(p- 1 ) T p (D) 



Since the operator U is not unitary, we do not use * but f for tp. The triangularity of the operator U 
guaranties that the Fock space H* constructed in terms of t/'iV' coincides with the original one. 
Thus, we can rewrite ( p7| ) as follows 

oo 
3 = 1 

The important property of this formula is that for given number of solitons n the summation can be taken 
from 1 to n because the operators ip2j-i with j > n produce "exact forms" when plugged into the matrix 
elements (see the proof of Proposition 2). 

Similarly, we can express the operator Q, defined in (Eq), in terms of ip: 



dD 
~D~' 



X(D) 



This equality is due to the fact that only the regular part of ip(D) contributes into the integral which does 
not change under the transformation to ip(D). 
Now we are ready to consider the operator Q>. 

Proposition 4'. From Proposition 4 one gets the following set of null-vectors due the vanishing of the 
residues: 

exp(^t 2fc _ 1 7 2fe _ 1 (B)) (*i|Q t exp(-^i/ ) ,_ 2fe J 2fc (A,B))| - 1} 



k>l 



k>l 



where 



Proof. 



Directly from Proposition 4 one gets the following formula for Q': 

Q f =res A=00 Y[^(A,B J )a- n J dD D 2n e~ y ^ V* \D) 



A 2 



(28) 



\D\>\A\ 



This formula looks much simpler in terms of xjr . By definition we have 



^*{D) = U T ^ (D) = 

\\Di\>\D\ 

The last term does not contribute to the residue because 



odd 



dDD 2n e~ Y ^ 



1 



D 2 -A 2 



^(D) smg = 0(A 2n - 2 ) 



\D\>\A\ 



and 



H*l>(A,B j )a- n =A- 2n (l + 0{A- 1 )) 



(29) 



Substituting the rest into ( j28| ) one has 

Qt = y^l^t (£)l)e X(^) reSA=oo ]^^ (A)jB . )a -n J dDP{D) _}__ T (R^ 



|Di|>|-D| 



13 



Using the formulae from Appendix B one can show that 



/ ^>tct (£) - 5 (Irk + i^k) + ° {A " ] = - 42 ""' 1 + 0(A ~'» 

|U 1 |>|£»|>|A| 

Here the equality is up to "exact form" in A; such "exact form" never contribute to the residue. Now the 
formula ( p9|) gives 

res A=00 J] ip(A, Bj )a~ n J d DP{D) ^-^ r e (j^ = 1 

j |Oll>l-D|>l^l 

which proves the proposition. ■ 



This alternative expression for Q< shows that it is independent of n, as Q and C are. 
Notice that in the formulae for Q and Q} only the holomorphic parts of ip(D) and tp^(D) are relevant. 
This leads to the important commutation relation: 

[C, 0)] = Q (30) 

Notice also that Q and are nilpotent operators, Q 2 = (Q^) 2 = 0, and [C, Q] = 0. 

This is the proper place to discuss the problems which we had before: the definition of local operator 
corresponding to arbitrary polynomial satisfying (^|) and the singularities in the definition of C and C e . 

Consider the polynomials 

- A 2 X*-3|Qexp(-£ ±h_ 2k J 2k (A,B))\ 1) (31) 

i<j k>l 

where the states (\&-3| are the vectors from the Fock space constructed via the fermions the ip. By the 
very definition of tp the polynomial (^lj) is equivalent to an antisymmetric polynomial ,{Ai, ■ ■ ■ ,A n ) 

of degree < 2n — 1 with respect to any Ai which is independent of the coupling constant. This polynomial 
satisfies the requirement formulated in and hence defines a local operator for arbitrary coupling constant. 
So, if we start from the fermions ip the counting of the local operators is independent of the coupling constant 
like in the paper [jjj . But where is the origin of vanishing denominators? The point is that if we consider the 
rational coupling constant £, for many of these local operators all the form factors will vanish. So, we have 
to be more careful: to consider the vicinity of the rational coupling constant £ = ir^ + e and to keep the first 
non-vanishing order in the form factors in the limit e — > 0. Thus, to define the space of local operators it is 
better to start from the fermions if). If when passing to the fermions if) we find somewhere infinite coefficient 
when the coupling constant is rational, it is always accompanied by a vanishing operator in such a way that 
the result is finite. 



3.3 Algebraic definition of the null-vectors. 

Let us summarize our study of the null- vectors for the descendents of <i>o. It is quite convenient to present 
the null- vectors as vectors from the dual space. In this section we shall write (^j = if the matrix element 
of (v^l with the fermionic generating function of local fields vanishes under the integral (w means "in a weak 
sense"). We have found three types of null- vectors: 

(*_ 5 |C =0, V (#_ 5 |GiP 5 , 

<*- 3 |Q =0, V (*_ 3 | eir 3 , 
(*i|Q f =0, V (*i|efff, 

where the operators C, Q\ Q are given respectively by 

C = J ^>U>U, Q = J ^-e^(D), Qt= J d -^e x ^{D) 
Let us show that these three conditions for the null- vectors are not independent. 



(ii) 
(Hi) 
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It is easy to show that the operator C identifies the spaces HZ 3 and iJj": 

Ker{C\ H *^ H ~) = 0, Im{C\ HU ^ H *) = H{ 

Hence every (^i | G HI can be presented as (\I/_ 3 |C for some (^-3! £ H*l 3 . Let us show that the null-vectors 
(iii) are linear combinations of (i) and (m). We have: 

(^|Qt = {^_ 3 |CQ t = {^- 3 \Q + (^-3\Q f C 

where we have used the commutation relation (pw. Thus we have proven the following 



Proposition 5. In the space of descendents of $0 = 1 which is H(I) (g> fP : (where -ff (I) is the space of 
polynomials of {hk-i}) the null-vectors coincide with the vectors 

(*_ 5 |C =0, V (*_ B | G#* 5 , (i) 
(*_ 3 |Q =0, V <*_ 3 | eHl 3 , (ii) 

and their descendents with respect to I 's. 

The consideration of the other operators $2m is based on the same formulae, but involves additional 
complications. We do not want to go into details, and we only present the final result. 

Proposition 6. For the operator $2m whose descendents are counted by the vectors from H(I) ® H^m-i 
we have two types of independent null-vectors: 

<*-5- 2m |(Cr +1 =0, V (*_ 5 _ 2ro | G #* 5 _2m, (i) 

(*-3- 2m |(C) m Q =0, V (*-3-3m| € #* 3 _2m, («) 

and their descendents with respect to I 's. 

For the operators <&2m+i one has the following picture. There is no uniform "exact form" which is an 
odd polynomial, so, the analogue of the operator Q does not exist for $2m+i, and the null- vectors are either 
due to the deformed Riemann bilinear identity or due to the vanishing of the residue. To construct the 
null- vectors in terms of fermions one has to introduce first the operator C: 

D 2 |>|-Di 

where the fermions are from the Ramond sector. This operators can be rewritten in a form similar to (|2 

C = / —■4){D) re g'4)(D) S i ng 

The fermions ip are related to tp by triangular transformation. 

The consideration of the operator $1 is absolutely parallel to the consideration of $0. The null- vectors 
are created either by the action of C (Riemann identity) or by the action of i/jq (residue) . Notice that 

[C,$]=0 

which guaranties the consistency. For higher operators $ 2 m+i there are additional problems which we would 
not like to discuss here. The general result is given in the following 

Proposition 7. The descendents of the operator $2m+i are counted by the vectors from the space H(I) ® 
H 2m . We have two types of null-vectors 

(*-4-2 m |(C) m+1 =0, V (*_ 4 -2m| G HU_ 2m , (i) 

<*2-2m|(C) m $ =0, V (* 2 _ 2m | G m_ 2m , (li) 

and their descendents with respect to I 's 
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3.4 Examples of null-vectors and the characters. 

Let us present the simplest examples of null- vectors for the operators $o, $1 and $ 2 - 
For the operator $ the simplest null-vector is created by 

(-3|G = *i(£)<-l| 

This null-vector is 

(iif)*" 

This null- vector is to be compared with 



L-i<Z>o 

For the operator $1 the simplest null- vector is created by 

-,2 1 



which gives the null-vector 



which has to be compared with 



(L_ 2 + kL 2 _ 1 )^ 1 

For the operator $2 the simplest null-vector is created by (—b\CQ. It yields 



which has to be compared with 



(L_ 3 + KiL_iL_ 2 + K 2 Lti)^ 
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Notice that the relative coefficients in our parametrization of the null-vectors are independent of £. So, they 
are exactly of the same form as the classical one. However, this is not always the case. 

Let us show that generally the number of our null-vectors is the same as for the representations of the 
Virasoro algebra. Recall that we consider the null- vectors which do not depend on the arithmetical properties 
of so, there is one basic null- vector in every Verma module of Virasoro algebra. The character of the 
irreducible module associated with <I> m is 

*'" <p) = (w '" +1, iwb*) (32) 

where we omitted the multiplier with the scaling dimension of the primary field. We can not control this 
scaling dimension, the dimensions of the descendents are understood relatively to the dimension of the 
primary field. The character ([32|) is obtained from the character of the Verma module by omitting the 
module of descendents of the null- vector on the level m + 1. 

Let us consider the character of the module which we constructed in terms of /, J. The dimensions of 
Iik-l a nd 3 2 k are naturally 2k — 1 and 2k. If we do not take into account the null- vectors, the characters of 
all the modules associated with <I> m are the same: 

x(p) 



IL>i(i-p*) 



Let us take into account the null-vectors. They are described in terms of fermions. By consistency with the 
dimensions of I^k-i an d J 2 k one find that the dimensions of ipi and ^Lz equal I. 

Technically it is easier to start with $ 2m+ i. The space of descendents is H(I) ® -ff| m where H{I) is the 
space of polynomials of {Izk-x}- 



1G 



Proposition 8. The character of the space of descendents of $2m+i, modulo the null vectors, equals 

X2m+ i(p) = (l-p 2(m+1) ) ' 



Proof. 

The null- vectors are defined in the Proposition 7. It is easy to eliminate the null- vectors (ii): we have to 
consider the subspace H*_ 2m o m which the level tpo is always occupied. Consider the sequence 

TT* C. tt* C7 1 rr* 

-"-2/71-4, ~^ n -2m, 11 2m, 

The operator C m identifies the spaces H*_ 2m and H 2m : 

Ker{C m \ HUm J = 0, Im{C m \ H ^ m J = H* 2m , 

and (—2m — 2|^>JC m = (2m|. Hence we can count the descendents by vectors of the space H(I) £g> H*_ 2m 
with the null- vectors: 

(*-4-2m|C =0, V (*_ 4 -2m|eir 4 _ 2mi , 

Notice that the operator C is dimensionless and 

Ker(C\ Hli _ 2mi ^i 2m , ) = 
That is why for the character of the space of descendents without null-vectors we have: 



X2m+1 (p) = / ( X H . (p) - X „. (P)) 



(33) 



where the first multiplier comes from H (I) the multiplier p m ( m ^ is needed in order to cancel the dimension 
of the vacuum vector in ff* 2m . Let us evaluate the expression in brackets: 

X H , (p) ~ X H * (p) = 

-2m, -2m-4, 

, dx 



x m - 2 - 



= f Y[(l+p 2j x)(l+p 2 ^x- 1 )x- m — - f Hil+p^x^l+p^x- 1 ) 
= J(l + x- 1 ) + p 2 ix)(l + p^x' x )x- m ^ — J (1 + x- 1 ) + p 2 ^)(l 
= {l-p 2{m+1) ) f (l + x- 1 )Y[(l+p 2] x)(l+p 23 x- 1 )x 

3 3>1 



27 -1\ -m^ X 

n J m 1 'T" — 



where we have changed the variable of integration x — > irp~ 2 in the second integral when passing from third 
to forth line. Substituting this result into (|33] ) we get the correct character: 

X2m+ i(p) = (l-p 2(m+1) ) ' 



Let us consider now the operators $2m- We parametrize the descendents of $ 2m by the vectors from 
Proposition 9. The character of the space of descendents of $ 2m , modulo the null-vectors, equals 

X2 m (p) = (l-P 2m+I ) fr ^ JT 

IL>i(i -P 3 ) 
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Proof. 

The null-vectors are defined in the Proposition 6. We have 

rr* tt* C™ TJ* 

TJ* Q TT* CT IT* 

rz -2m-3 -"-2171-1 fl 2m-l 

The operator C m identifies H*_ 2m _ x and i?|m-i : 

and {—2m — 1 C m = (2m — 1|. Hence we can replace the space H{I) ® ff|m-i these null-vectors by 
_ff (I) <g> ff* 2m _i with null-vectors 

(* -5-2m|C = 0, V (*_ 5 -2m| € #* 5 _ 2m , (*) 

(*-3-2m|C=0, V (*_ 3 -2m| € -ff-3- 2m! (**) 

So, the character in question is 

X2m{p) = Yj 7^ 27TTTP~ m2 (x „. (?) - X H . (p)) 

llj>ivl _ P J ^ -2m-l, -2m-5,0 / 

where -ff* 2 ;-i o = ^-21-1/^-21-3^- ^ n order to calculate the character x . (p) one has to take into 

' -21-1, 

account that Q is a nilpotent operator, Q 2 = , with a trivial cohomology. Hence 

Mfikv,^.^) - MQk- a ,_ B ^- a ,_,) 

Summing up over this complex we obtain: 

dx 

i+p" J *x)(l +p" J -^x-^)x~" 

T. -4- 1 

J>1 



x,. (p)= / nti+p^ 11 )^^ 1 ^ 1 )^^- 

-27-1, J - t J - X + 1 X 



|X|>1 

Hence 



X H . (p) X jj* (p) = 

-2m-l, -2m-5, 

J -^1 .T + 1 



|x|>l ^ 



f + p 2 i- l x)(l + p 2 i- l x- l )x- m — - f \{(l+p 2 i- l x)(l+p 2 i- l x- l )x- m - 1 — 



J>1 " J>1 

= (i-P 2m+1 )x^ (p)=P m2 (i-P 2ro+1 ) 1 



n,->i(i-p 2j ) 

Thus the character is given by 

X 2m (p) = (l-p 2m+1 ) 1 



n J >i(i-p j ) 

as it should be. 



4 Classical case. 

4.1 Local fields and null- vectors in the classical theory. 

The classical limit of the light-cone component T__ of the energy-momentum tensor gives the KdV field 
u(x-). When working with the multi-time formalism we shall identify X- with t\. Local fields in the KdV 
theory, descendents of the identity operators, are simply polynomials in u(t) and its derivatives with respect 

to t\\ 

O = 0{u,u',u",...) (34) 
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We shall use both notations d\ and ' for the derivatives with respect to X- =t\. 

Instead of the variables u, u', u" , it will be more convenient to replace the odd derivatives of u by the 
higher time derivatives d 2 k-\u, according to the equations of motion 



0L 
dt 2 k-i 



2fc-l 

L— 



,L 



1 



2 2fe-l 



,(2fc-l) 



Here 

L = d\ - u 

is the Lax operator of KdV. We have used the pseudo-differential operator formalism. We follow the book 
The even derivatives of u(x) will be replaced by the densities S^fe of the local integrals of motion, 



S 2k = res dl L^ = -——. u (^ 2 ) + • • • , 



2 2fe- 

In particular on level 2 we have S 2 = ~\ u - For a reader who prefers the r-function language S 2 k = 
did 2k -i logr. 

From analogy with the conformal case we put forward the following main conjecture underlying the 
classical picture: 

Conjecture. We can write any local fields as 

0(«,«', «",...) = F O:0 (S 2 ,S 4 ,---) + J2d u F o ^(S 2 ,S 4 ,---) 



(35) 



where v = i 3 , ■ ■ ■) is a multi index, d v — 1 1 ; 



We have checked this conjecture up to very high levels. To see that this conjecture is a non trivial one, 
let us compute the character of the space of local fields eq. (p4|) . Attributing the degree 2 to u and 1 to 0i, 
we find that 



n t~~~i = c 1 -p) n T~~~i = 1 +p 2 +p 3 + 2 p 4 + 2 p 



J>2 r j>l 

On the other hand the character of the elements in the right hand side of eq. ((3^) is 

X2 



n r^rr n = n d^j = 1 + p + 2p2 + 3p3 + 5p4 



7 P 5 



j>i i>i j>i 

Hence xi < %2, and this is precisely why null- vectors exist. Let us give some examples of null- vectors 



level 1 


01 • 1 = 










level 2 


0? ■ 1 = 










level 3 


0? ' 1 = 0, 


03 • 1 = 








level 4 


0* • 1 = 0, 


0103 • 1 = 


0, 


(dfS 2 - 45 4 


+ 6Sf ) -1 = 


level 5 


0i 5 • 1 = 0, 


0! 2 03 • 1 = 


o, 


05 • 1 = 0, 






01(0^2- 


4^4 + QSl) 


1 = 


0, (d 3 S 2 " 


0154) -1 = 



(36) 



We have written all the null-vectors explicitly to show that their numbers exactly match the character 
formulae. The non trivial null- vector at level 4 expresses S4 in terms of the original variable u: 4S4 — 
— \u" + |w 2 . With this identification the non-trivial null- vector at level 5, d 3 S 2 — 0i<?4, gives the KdV 
equation itself 

3 u + ~mt' - hul" = 

More generally one can consider the descendents of the fields e miv where ip is related to u by the Miura 
transformation 

u = —ip + up 
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Here, the presence of i is a matter of convention. The reality problems have been discussed at length in 

For this consideration and for other purposes we need certain information about the Baker- Akhiczcr 
function. The Baker- Akhiezer function w(t, A) is a solution of the equation 

Lw(t, A) = A 2 w(t, A) (37) 

which admits an asymptotic expansion at A = oo of the form 

w(t, A) = e^' A \l + 0(1/ A)); C(*, A) = £ h^A 2 ^ 1 

k>l 

In these formulae, higher times are considered as parameters. The second solution of equation (|37]), denoted 
by w*(t,A), has the asymptotics 

w*(t,A) = e _c(t ' A) (1 + 0(1/4)) 

These definitions do not fix completely the Baker-Akhiezer functions since we can still multiply them by 
constant asymptotic series of the form 1 + 0(1/ A). Since normalizations will be important to us, let us give 
a more precise definition. We first introduce the dressing operator 



i>l 

and we define 

w(t, A) = Se^, w*(t, A) = ($*)- 1 e -«*'^ 

where $* = 1 + X),>i( — di) ~*<I>i is the formal adjoint of 
Proposition 10. With the above definitions, one has 

1 ) The wronskian W(A) = w(t, A)'w*(t, A) - w*(t, A)'w(t, A) takes the value 

W(A) = 2A 

2) The generating function of the local densities S(A) = 1 + J2k>o ^2hA~ is related to the Baker-Akhiezer 
function by 

S(A) = w(t,A)w*(t,A) 

3) The function S(A) satishes the Ricatti equation 

2S(A)S(A)" - (S(A)') 2 - AuS(A) 2 - AA 2 S(A) 2 + AA 2 = (38) 

Proof. 

Let us prove the wronskian identity. This amounts to showing that resA(W(A)A l ) = 2<5, l _2. But we have 

res A (w(A)A^ = res A { (ftSd'e^^) (V)-^ 4 ^) - {^e^ A Ah x (^*)-\-d 1 fe-^ t ' A ^\ 
We can transform the residue in A in a residue in d\ using the formula 

res A { (Pe«^)) • (Qe~^) } = res dl (PQ*) 

Hence we find 

res A (w(A)A i ) = resg 1 |9i<i>9j$ _1 + <f>9j$ _1 3i| = resg 1 ^d\L^ + L^d\^ 

If i is even positive, the residue is zero because the is a purely differential operator. If i = —2 the residue 
is obviously 2, and if i < —2 it is zero. If i is odd, then (La)* = -Ls so that the operator d\L^ + L^d\ is 
formally self-adjoint and it cannot have a residue. 
The proof of 2) is simple p3|| 



res a 



(A 2k - 1 w(t,A)w*(t,A)) =res dl ($df k - 1 $- 1 ) = res 9l (l,^) = S 2k 



20 



The Ricatti equation follows immediately from 1),2) and eq.(|37|). 



Let us return to the descendents of the primary fields. For the descendents of the fields e mv> , our 
conjecture states that 



(39) 



u>0 



Let us consider for example e zip . For a true solution of the KdV equation, the Baker- Akhiezer function is 
a true function on the spectral curve, and it can be analytically continued at A = 0. From the definition of 
e lip we have Le l(p = 0. Comparing with eq.(^), we see that e l(p — w(t, A)\a=q- To check eq.(^), at least on 
the first few levels, we need the time derivatives of e llp . They can be obtained as follows. The time evolution 
of the Baker- Akhiezer function is well known. 



dw 
dt-2k-i 



w 



By analytical continuation at A = 0, we obtain the evolution equations for e ltp . 

Let us give some examples of these null vectors. We show below the first null vector associated to the 
primary fields e lmv> . 



m = 1 

m = 2 

m = 3 

m = 4 



(df + 2s 2 y v = o 

(2d 3 + Of + &d x S 2 )e 2iv = 

(8<9i<9 3 + df + 12d 2 S 2 + 24S 4 )e 3itp = 

(24<9 5 <9i + 2093d? + df + 20<9 1 S' 2 + 40939152 + 120dfS 4 )e illp = 



In these formulae, the derivatives act on everything on their right i.e. diS^e ^ = d\(S2e 



2i V \ 



4.2 Finite-zone and soliton solutions. 

For the finite-zone solutions, the Baker- Akhiezer function is an analytical function on the spectral curve 



which is an algebraic Riemann surface. Let us recall briefly the construction 14, [L5| 



We start with an hyperelliptic curve T of genus n described by the equation 

2n 

T :Y 2 = XV(X), P(X) = l[(X-B 2 ), B 2n > •••> B 2 > B!>0 

For historical reasons we prefer to work with the parameter A such that X — A 2 . The surface is realized 
as the A-plane with cuts on the real axis over the intervals c; = (B 2 i~i, B 2 i) and q = (—B 2 i,—B 2 i-i), 
i = 1, • • ■ , n, the upper (lower) bank of Cj is identified with the upper (lower) bank of Ci. The square root 
yJV{A 2 ) is chosen so that y/V(A 2 ) — > A 2n as A — > cxo. The canonical basis of cycles is chosen as follows: the 
cycle di starts from i?2i-i and goes in the upper half-plane to — i?2«-i, &i is an anti-clockwise cycle around 
the cut Cj. 

Let us consider in addition a divisor of order n on the surface T: 



V=(P 1 ,---,P n ) 

With these data we construct the Baker- Akhiezer function which is the unique function with the following 
analytical properties: 

• It has an essential sing ularity at infinity: w(t, A) = e c(t ' A) (l + 0(1/ A)). 

• It has n simple poles outside infinity. The divisor of these poles is T>. 

Considering the quantity —d 2 w + A 2 w, we see that it has the same analytical properties as w itself, apart 
for the first normalization condition. Hence, because w is unique, there exists a function u(t) such that 

- dfw + u(t)w + A 2 w = (40) 
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We recognize eq.(|37]). One can give various explicit constructions of the Baker- Akhiezer function. Let us 
introduce the divisor Z(t) of the zeroes of the Baker-Akhiezer function. It is of degree n: 

Z(t) = (A 1 (t),---,A n (t)) 
The equations of motion for the divisor Z(t) read |D| . 

dA (t)- v^W) (ad 

diMt) -umrm) (41) 

The normalization of the Baker-Akhiezer function corresponds to a particular choice of the divisor of its 
poles V. Later we shall specify the divisor which corresponds to the normalization of the Baker-Akhiezer 
function which was required in the previous subsection, for the moment we give a formula in which the 
normalization is irrelevant. Consider two sets of times t and i' ', differing only by the value of t\. Then we 
can write 



W(t ' A) -.l Q W> exp r^S^x (42) 



w{m,A) y Q(A 2 ,m) \J t w Q(A 2 ,t) 

where the polynomial Q(A 2 ,t) is defined as Q(A 2 ,t) = JJ^A 2 - A 2 (t)). 

The ratio of two dual Baker-Akhiezer functions ^t^ta) * s obtained by applying the hyperelliptic invo- 
lution. This amounts to the reflection A — > —A in eq. (|42|) . Let us prove the following simple proposition. 
Proposition 11. For the Baker-Akhiezer functions w(t, A), w*(t, A) normalized by 

w(t, A)'w* (t, A) - w* (t, A)'w(t, A) = 2A 

we have 

5(A) = -^j= ^ OX], I > "f/.^.l - I (4.T) 



the latter equality is the definition of J^k- We recall that Q(A 2 ) and V(A 2 ) are the polynomials 

n 2n 

Q(A 2 ) = Y[(A 2 ~A 2 ), T(A 2 ) = l[(A 2 -B 2 ) 

i=l z=l 

Proof. 

To prove the proposition we use the Wronskian identity 

w(t,A)'w*(t,A) -w*(t,A)'w(t,A) = w(t,A)w*(t,A)d 1 ( log W ^ ,A }, '] = 2A 

\ w*(t,A)J 

but using eq.(f42|) and the fact that ^/(o'^l) and ^j0Ta) differ by the sign of the square root, we have 



w*(t,A)J Q(A 2 ) 



and the result follows. 



Notice that for J 2 k defined in ( p3| ) we have 

i i 

From eq.([43|) we see that the normalization of w(t, A) and w*(t,A) which corresponds to the proper 
value of the Wronskian is such that the divisors T> and V* are composed of Weierstrass points and V + V* = 
{Bi, ■ ■ ■ , B 2n ). Actually, it is this quite unique normalization which was used by Akhiezer in his original 
paper. 

Now we are in position to describe the dynamics of S(A) with respect to all times. It is very useful to 
define the following strange object 

dl(D) = V D- 2k —^—dD (44) 

fe>l Zh 1 

dI{D) is a 1-form in the Z?-plane and a vector field with respect to times. We have 
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Proposition 12. 



turn- s ( A) g ww-yw iD (45) 

Proof. 

We give a proof of this proposition for the finite zone solutions, which are our main concern here, but clearly 
the formula is quite general. We are sure that a general proof of equation ([15]) exists, but it must be based 
on manipulations with asymptotic formulae. We prefer to work with analytical functions. Anyway, every 
solution of KdV can be obtained from the finite-zone ones by a suitable limiting procedure, so considering 
finite-zone solutions is not a real restriction. 

Let us describe the motion, under the time ij, of the divisor Z(t) of the zeroes of the Baker- Akhiezer 
function. Introduce the normalized holomorphic differentials du>i for i = 1, • • ■ ,n and the normalized second 
kind differentials with singularity at infinity du2%—i, i > 1 

dwi = Si j 

J du 2l -i = 0, dw 2l -i(A) = d{A 21 - 1 ) + 0(A- 2 )dA for A ~ oo 
It is well known that, by the Abel map, this motion is transformed into a linear flow on the Jacobi variety. 



A* 



^ / duj k = I dio 



2Z-1 



dt 2 i-i 

3 

Form this equation one easily find : 

dI{D) ■ ( ' dcj k = ( dZ) D (46) 

where duju is a 2-differential defined on T x T (r is the Riemann sphere) parametrized by A and D 
respectively. It is useful to think of T as a realization of the curve Y 2 = X similar to T. The a-periods of 
duii) on r vanish. The only singularities of the differential duiu are the second order poles at the two points 
A = ±D: 

d ^{A) = ( ,f a *fl 8)a + 0(X))dAdD 
By Riemann's bilinear relations one easily find: 

/ dui D = du k (D) 

Eq.(^) then takes the form 

dl(D) -Y, [ 3 dcok= dcu k (D) (47) 

3 

The normalized holomorphic differentials arc linear combinations of the differentials 

da k (A) = - dA, k = l---n 

with coefficients depending on Bi. They do not depend on times. Hence by linearity we can write for da k 
the same equation as (|47|). Differentiating explicitly we get the following system of equations 

n j^2k — 2 jy2k—1 

V 4 dim) ■ A, = dD; k = 1 • • • n 
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Solving this linear system of equations gives 



Q(D 2 ) i y/Hty Jn _ Q(D*) I 



dI(D) ■ A, = 2A j ^ y ' — -j= v . A0 -r^-dD = — ^ v ' — -rd x A i dD 

3 3 v^W) D ~ A \ 'iW^-4 ? ) jnp^Di-A) * 

where we have used eq.(p[) in the last step. Finally, using eq.([43|) we find 

dl(D) ■ S(A) = S(D)S(A)^2 D2 l _ A 2 A 2 l _ A M A])dD = ^-_S(D)S(A) (l°g§^)' dD 



The soliton solutions correspond to a rational degeneration of the finite-zone solutions such that 

B 2 j-i — > Bj <— B 2 j, j = 1, • • • , n 

The points of the divisor At and the points Bj are the coordinates in the n-solitons phase space. In wc 
gave a detailed discussion of the Hamiltonian structure. In particular we have 



- = {hk-i, •}; hk-i - 2k _ 1 ^2 B f k 1 



dt 2k . 

The expressions for the local observables also follow easily from the finite-zone case. In particular 

n 

J 2fc (A,5)=£(Af -Bf) (48) 

i=l 

The expressions for S 2 k follow from here. 
4.3 Classical limit of Q and C. 

Let us consider the classical limit v — > oo of the operators Q and C. To this aim, we have to understand 
the relation between the quantum and the classical descriptions of the observables. In the quantum case we 
considered the form factors i.e. matrix elements of the form 

fo(/3i, ■ • ■ , &»)_..._+...+ = (010(0)1/?!, ••■,/?„; n+1 , • ■ • , f3 2n ) 

where •■•,/?„ are rapidities of anti-solitons, n +i, ■ • • , Pin are rapidities of solitons. The matrix ele- 
ments of this form do not allow a direct semi-classical interpretation, it is necessary to perform a crossing 
transformation to the matrix elements between two n-soliton states: 



(Pi, ■ ■ ■ , n \O(O)\j3 n+ i, ■ ■ -,P 2n ) = foWi ~ m, ■ ■ ■ , Pn ~ ni,/3 n +i, • • -,#2n)- 



__i — h 



In |Q it is explained that the formula (Q) for this form factor is a result of quantization of n-soliton solutions in 
which Ax, ■ ■ ■ , A n play the role of coordinates, B\, ■ ■ ■ , B n and B n +\. ■ • ■ , B 2n give the collection of eigenvalues 
for two eigenstates. 

Recall that the generating function for the local descendents of the primary field $ m was written as 
follows 

(n 2n \ 

n^n^ 7 * (49) 
i=l 3=1 ) 

where I 2 k-i{B) and J 2 k(A\B) arc defined in eqs.(|]j|). The expression 



/ n 2n 

exp(£ y 2k J 2k (A\B)) ]J AT ]J B j 



k>l \i=l j=l 
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is practically unchanged under the crossing transformation which corresponds to Bi — * —Bi for i = 1, • • • , n. 
Comparing it with the classical formulae ( f48| ) we see that it corresponds to special symmetric ordering of 
them, for example 



J2k{A\B) — 2 (j2k(A%, ■ ■ • , A n , Bi, • ■ • , B n ) + J2 k {Ai, • ■ ■ , A n , B n+ i, 



Bo 



This ordering is a prescription which we make for the quantization. 

On the other hand the eigenvalues of the Hamiltanians l2k-i{B) under crossing transformation change 

to 

—l2k-i(Bi, • • • , B n ) + l2k-i(B n+ i, • • • , B 2n ) 

i.e. the descendents with respect to /2fc-i correspond to taking commutator of O with /2fc— i- Certainly the 
classical limit makes sense only for the states with close eigenvalues, so, it is needed that 

S2k-l{Bl, • • • , B n ) — S2k-l{B n +l, ■ ■ ■ , B2n) = C(£) 

Recall that £ = — ilog(g) plays the role of Plank's constant. 

Thus comparing the classical and quantum pictures provides the following result. The quantum gener- 
ating function (^9|) corresponds to the classical generating function 

£*(t, y) = exp(j2 *a*-i Ja*-i) • exp(£ y 2 k J 2fc ) (50) 
fe>i fc>i 

where • means the application of Poisson brackets. In fact /2/t-i can be replaced by &2k-i- The normalization 
in the formula for 1(B) (^) is chosen in order to provide an exact agreement with the classical formulae. 

Now let us consider the classical limit of the operators Q and C in the Neveu-Schwarz sector. For Q we 
had the formula 



Q = J ^^)f(fl) = / ^smh(X(DmD) 



The latter equation is due to the fact that the fermion is odd. From the definition (22) of X(D), we have in 
the classical limit 



k>l 



'2k-l 



Hence the following expression is finite in the classical limit: 

Qa = lim i2j = J yj(D)dI(D) (51) 
where dI(D) is the 1-form in Z?-plane introduced in the previous subsection: dI(D) = J2 D- 2k I 2k ~idD. 

k>l 

Remark that Q c i can be thought of as a generalized Dirac operator. 
For C we had the formula 



dDj 

£>2|>|-Dl| 



' Uh 1 " lh cosh(Xp 1 ))cosh(X( J D 2 ))r- I ^{D^{p i} 



C = 



D 2 J Di v v " v v " e \D 2 

D 2 |>|-Di| 

J ^smh(Xpx))smh(X(D 2 ))r e + (g) V>0W(AO 

|D 2 |>|-Di| 

where t+ and t~ are even and odd parts of r e : 

,2fc-l 00 1 , „2fc 



W = 2^ i t 2 fe-i x > w) = -Lrr^ 



x 2k 
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Obviously when £ — > one has 



ret*) - ( -— ) , T+( X ) -> -(i^ log(l - Z 2 ) 



2 da; VI 

So, the following expression is finite in the classical limit 
Cci = lira — = 



J \p(D)^(D)dD + | dI{D 2 ) J dialog U 

^D 2 \>\D^ 



,02 



il>(Dx)il>(D 2 ) (52) 



In the next subsection we are going to apply these operators to description of the classical KdV hierarchy. 
Notice that as usual the quantum formulae are far more symmetric than the classical ones. 



4.4 The classical equations of motion from Q d and C, 



ch 



In this subsection we shall consider only the descendents of the identity, i.e. the pure KdV fields. We have 
described this space by the generating function (E(J): 



£m=o(*>2/) = exp^^i 2 /c-i^2fe-i) -exp(^y 2 fc</2/c) ■ 1, 



k>l fe>l 
Let us fermionize J 2 k and apply the equations 

(*- 3 |eoi = 0, (*- 5 |C ci = (53) 

to the description of the equations of motion. In this section the symbol = means the vanishing of the 
scalar product with the generating of local fields. 

We give the list of null-vectors following from these two equations up to the level 5 specifying explicitly 
the vectors (^-3) and (vfr-sl from which they come. We do not write the descendents with respect to J's of 
the already listed null vectors: 
Null vectors coming from Q c ;: 



d 5 -l 

(-d 3 S 2 + 9i5 4 ) • 1 



(-IIVI3C1 : (dfS 2 - 4S 4 + 6S% + \d x d 3 ) • 1 



Null vectors coming from C c 



Obviously these null- vectors coincide with (p6[). So, in particular, equations ( |53| ) imply the KdV equa- 
tion itself. We have verified that the null-vectors coincide with those obtained from the Gelfand-Dickey 
construction up to level 16. On higher levels we find higher equations of KdV hierarchy. 

We have seen that the KdV equation follows from the equations (|5^). Let us prove the opposite: equations 
(|53]) hold on any solution of KdV. We start with the operator Q c i. 



Proposition 13. Let 



Qd = 



[ iP(D)dI(D) = 



k>l 

Then if J 2k are constructed from a solution of KdV we have 



d 

ohk- 



Qd exp - > - J 2k h-2k 1-1) = 
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Proof. 

Let us introduce the notation: 



T = exp (- g \j 2k h^ = exp (J 



logS(A)/»(A) 



where h(A) = £ h~ 2k A 2k . We have i/>(D) T = T S'^D^D). So 

k>l 

We now use eq.(|45|) to get 



/ 



i ^4 

-(dI(D)-S(A))\-l) 



S(D)S(A) Ty ' y '2iirA 
dA 1 
2inA D 2 - A 2 



HA 1 

dD ° n2 [(logS(A))' - (logS(2?))'] - 1) 



But 



Let us consider first the integral 

) . T AD 1 

1) 



dD dA D n 

(logS (A)) 



D{>{Al D 2iirAD 2 - A 2 



A D 

i,(D)i,(A)r(A):- w —^il,{A) 



One can do the integral over D. Notice that the integrand is regular at D = 0. Hence, the contributions 
to the integral come from the poles at D 2 — A 2 . The simple pole does not contribute because its residue 
vanishes since we have the product of two fermion fields at the same point in the normal product. The 
double pole obviously does not contribute either, and the integral is zero. Next we look at the integral 



dD dA D , 
(log S(D)) 



{D>lAl D 2mAD 2 -A 2 



A D 



1) 



This time one can do the integral over A. But it is clear that the integrand is regular at A = 0, and the 
integral also vanishes. ■ 



Let us consider now the operator C c i ■ 
Proposition 14. Let 

Cd = / dDyj(D)—ij(D) + 2 



'\Dx\<\ 

then if J 2k are constructed from a solution to KdV 



- [ log (l - ^) ^j(D 1 )^(D 2 )dI(D 1 )dI(D 2 ) 



C c/ exp ( - \ j 2 kh- 2k I 



-1) = 



Proof. 

Let us split C c i into two pieces C c i — C\ + C 2 : 



- / log (l - ^| ) it{D 1 )^(D 2 )dI{D 1 )dI{D 2 ) 

m J\D 1 \<\D^\ V D 2 ) 



° 2 2m 
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We use the same notation T as in the previous proposition. We treat first the C 2 term 



C 2 T | - 1) = ~J dI{D 2 )^{D 2 ) T J 



dD x dAi 
Di 2inA 1 



log 



5(A!)V 



|Ai|<|Ui|<|D 2 | 



£) 2 V £>?-A? V & S(D 1 ) 
AxD x 



D\-A\ 



1-1) 



by the same argument as in the previous proof, we see that only the term (log S(Ai))' contributes. This 
time however, the double pole gives a non vanishing contribution 



1 

2ni 



dD\ log 1 



|Z?i|>|Ai 



D\) {D\-A\Y ~ D\-A\ 



so that 



Co T I - 1) 



dI{D 2 )iP(D 2 ) T 



\D 2 \>\Ai\ 



2^ £>| — Af ^(^)l- 1 ) 



Recall that in this formula dI(D 2 ) acts on T and on (log S(Ai))'. Commuting again ip(D 2 ) and T we finally 
get after having computed the integral over D 2 



i(logS(A))" + i((logS(A))') 2 



The calculation of the C\ term is straightforward 

C X T\-1) = | j dAS-^A)^)-^^-!) 

Combining the two terms and using the Ricatti equation ( |38| ) 

2(log S{A))" + ((log 5(A))') 2 + AA 2 S- 2 {A) = 4(« + A 2 ) 

we get 

C cl T\-l) = j ^(u + A^(A)^(A)\-l) = 
the latter integral vanishes since the integrand is regular. 

The equation ( |53| ) for the generating function £ c '(i, y), defined by 

£ c '(t, y) = exp j — ^2 V2k-hk{t) 

can be written explicitely by bosonizing the fcrmions. The equation ("I" — 3 1 = then reads: 

j D e« D »*) dI(D)C cl (t, y - [D\ e ) = 
where dl(D) = J2k D~ 2k d 2 k-\dD as before (it acts on C cl by cfefc-i) and, 



fe>i 



The (\I>_5|C C ; = equation reads: 

87ri /" dD# 3 e 2i{D ' y) C c \t,y~2[D] c 



(54) 
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+ J D^Dl-DDlogU-^j e^ D ^+^ D ^dI(D 1 )dI(D 2 )C c \t,y-[D 1 ] e -[D 2 ] e )=0 
|£>a|>|Z»i| 

Equations ( |54| , |55| ) provide a system of linear differential equations for the Taylor coefficients of C cl (t,y). 
It becomes a system of non-linear differential equations for the J 2 k when the closure condition C cl (t,y) = 

exp fefe>i U2kJ2k(t)^j is imposed. 

Thus we have shown that the equations ( |53| ) provide the complete description of the KdV hierarchy. We 
find it quite amazing. It is also interesting that this new description came from pure quantum considerations. 



5 Connection with the Whitham method. 

There is an surprising relation between the methods of this paper and the Whitham equations for KdV 
[l7| , [Ts| , [nj . The present section is devoted to the description of this relation. 

Let us remind briefly what is the Whitham method about. Suppose we consider the solutions of KdV 
which are close to a given quasi-periodic solution. The latter is defined by the set of ends of zones , • • • , B 2n . 
We know that for the finite-zone solution the dynamics is linearized by the Abel transformation to the Jacobi 
variety of the hyper-elliptic surface Y 2 = XV (X) for V(X) = ]J(X - The idea of the Whitham method 
is to average over the fast motion over the Jacobi variety and to introduce "slow times" Tj which are related 
to the original KdV times as Tj — etj (e <C 1), assuming that the ends of zones Bj become functions of these 
"slow times" (recall that the ends of zones were the integrals of motion for the pure finite-zone solutions). 

For the given finite-zone solution the observables can be written in terms of ^-functions on the Jacobi 
variety, but this kind of formulae is inefficient for writing the averages. One has to undo the Abel trans- 
formation, and to write the observables in terms of the divisor Z = {A\, ■ ■ ■ ,A n ). The formulae for the 
observables are much more simple in these variables, and the averages can be written as abelian integrals, 
the Jacobian due to the Abel transformation is easy to calculate. The result of this calculation is as follows 
p7| . Every observable O can be written as an even symmetric function Lq(Ai, • ■ ■ , A n ) (depending on B's 
as on parameters), and for the average we have 



where 



A = det 



A^~ 2 dA 



1,3- 



The similarity of this formula with the formula for the form factors ([j]) is the first intriguing fact. Indeed, 
we have the following dictionary. For the local observables, we have 



For the weight of integration, we have 

I 



3=1 



But the most striking feature is that the cycles of integration are replaced by functions of a% = A\ v 

cycle ai a^ 1 

The coincidence between the notations for ai- variables and a^-cycles is therefore not fortuitous. The expla- 
nation of the fact that the cycles are replaced by these functions is given in [0] , where it was shown that the 
factor Yli <\ l selects the classical trajectory in the semi-classical approximation of eq.(^|). So, the solution of 
a non trivial, full fledged, quantum field theory has provided us with a very subtle definition of a quantum 
Riemann surface. 
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The main result of the Whitham theory is that the averaged equations of motion can be written in the 
following form: 

— du 2q+ i(A) = — du 2p+ i(A) 

<J±2p+l Ol2q+l 

where, as earlier, duj2 q +i, du2p+i are normalized second kind differentials with prescribed singularity at 
infinity: 

du 2l +i{A) = d(A 2l+1 ) + 0{A- 2 )dA, J oED 2/+1 = 

In fact every equation of this type contains many partial differential equations for Bj because one can 
decompose it with respect to the parameter A. To our knowledge these equations have never been deduced 
directly from the averaging integrals except for the special case p — 1 , q — 2 which was considered in the 
pioneering paper |17|] . We want to show that the methods of this paper allow to do that, also we shall also 
explain how to describe other equations of the hierarchy. 

As it has been shown in this paper the classical observables can be obtained from the generating function: 

exp(J^ tzk-ihk-i) ■ exp(J^ y 2fe hk) (55) 

k>l k>l 

In the Whitham case the KdV times become fast which means that l2 P -i = e^2p-i where / 2p _i denotes 
the Poisson bracket with corresponding Hamiltonian of the averaged hierarchy. To keep the first multiplier 
in (|5^) finite one has to pass to the slow times T2fc-i- So, the observables for the averaged hierarchy are 
generated by 

{(C (T,y\B))) =exp(^T 2fc _ 1 / 2fe _ 1 ) • (( exp(^ y 2k J 2 k) » 

k>l k>l 

Let us repeat that the average in this formulae depends only on Bj which are supposed to be functions of 
slow times, the multiplier with l2k-\ corresponds to derivations with respect to the slow times. 
The averaged formulae becomes really beautiful with the fermionic generating function: 

«£ (S)))| - 1} = A- 1 g(B) • • ■ «^»n| - 1 - 2n) (56) 

where 

g(B) is the same as in quantum case (fL7|). 

Remind that we have rewritten the equations of KdV hierarchy in the weak sense as follows 

<*_ 5 |C ci -0, (*_ 3 |Q C , = 

For the averaged hierarchy these equations have to be replaced by 

<*_ 5 |C -0, (*- 3 |Qo = (57) 

In this section, the symbol = means vanishing of the matrix elements with the averaged generating function 
(|56l ) . The averaged Co and Qo are 



Co = J 4iD)^{D)dD, Q = J 4>{D)dI{D) 



where the terms with /(D) are omitted in the definition of Co comparing with C c ; ( p2\ j because dl(D) is of 
order e, the definition of Qo is practically the same as in ( |5l|) because the latter is homogeneous in dl(D). 
Let us explain the implications of the equations (^) for the averaged hierarchy. 

Proposition 15. The equation 

(*_ 5 |C = 
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for the averaged hierarchy follows from the Riemann bilinear relation for hyper-elliptic integrals. 
Proof. 

The Riemann bilinear relation for the hyper-elliptic integrals is equivalent to the formula 

f dAi r dAa 

Cl c 2 

where c\ o c 2 is the intersection number of cycles, the anti-symmetric polynomial C(A±, A 2 ) is given by 

C(A 1 ,A 2 ) = C'(A 1 ,A 2 )-C'(A 2 ,A 1 ) 

where 

(see, for example, JO] for a relevant discussion). 
On the other hand we have 

Co g(B) = g{B) C 

where 

C = J P(D 2 )D- in ^(D)^(D)dD 

Consider the formula for averaged observables (|56|). Suppose that we undo the averaging considering instead 
of 

<*_ 1 |(^*» 1 .--(^*» n |-2n-l) 

the polynomial 

• • • r(A n )\ - 2n- 1) [] 

The averaging corresponds to integrating over closed cycles on the surface, so the latter polynomial is defined 
up to exact forms of the following kind 

f>ir m{a\, . ■ ■ , ai . ■ ■ , ad JhaJ)^ Uvwwm)) 

where Q(A) is an arbitrary polynomial and M(Af, A^ l _ 1 ) is anti-symmetric. 

It is clear that the statement of the proposition will be proven if we show that for every (^-5! up to 
exact forms one has 

(^. 5 \C r(Ax) ■ ■ ■ V{A n )\ - 1 - 2n> n A r _1 - 

j 

* ^(-l)^M(A?, ADCiA^Aj) (58) 

i<j 

for some anti-symmetric M{A\, ■ • • , A\_^). 

Just like in the proof of Proposition 2 from Section 3 we have three possibilities for (^-sl which can give 
a non-trivial result: 

1. The depth of is greater than —2n — 1 

2. The vector (^-sl is obtained from a vector (^-ll whose depth is greater than — 2ri — 1 by application 
of 2~ 1 V'-2g— 1 with q > p > n (i.e. there are two holes below — 2n — 1). 

3. The vector (^-sl is obtained from a vector (^-3! whose depth is greater than — 2n — 1 by application 
of ^* 2 p_i with p > n (i.e. there is one hole below — 2n — 1). 

In the first case using the identity 

r ^/vm d 

J D 2 -A\ dD 

|0|>|Ai|,|.A 3 | 
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one easily gets the formula ( ]58| ) with 

M(Al ■ ■ -,A 2 n _ 2 ) = (*- 5 \r(Ai) ■ ■■r(A n - 2 )\ 1 - 2n> JjAf " x 

In the second case it is necessary that in the expression (^-il^* 2 j,_iV ; *-2g-i the operator Co annihilates 
two holes. Hence one find the integral 



(*_i|2(p - g) /" V(D 2 )D~ An+2p+2q+1 dD = 



It vanishes because p,q > n. 

Finally, in the third case it is necessary that in the expression ^s\ip^_ 2p —i Co the operator Cq annihilates 
the hole. This gives 

<*_ 3 | J (v(D 2 ) + D^P(D 2 )^j D- in+2 ^(D)dD 
So, in the matrix elements we shall find the polynomials 

\D\>\Aj\ 3 

which corresponds to an exact form. ■ 

The similarity of the proof of this proposition with that of Proposition 2 of Section 3 is quite impressive. 
Let us consider now the operator Qq. It is responsible for the equation of motion as shows the following 
proposition. 

Proposition 16. The equations 

d d 
— du 2q+l (A) = — <£j 2p+1 {A) (59) 

Ol2p+l U±2q+1 

follow from 

(*- 3 |Qo = <*_ 5 |C = (60) 

Proof. 

We will show that the Whitham equations (^9|) follow by considering the vectors 

(*_ 3 | = (-11^*2^-1^*29-1^+1 

The proof goes in two steps. First we shall show that the equation (|6(]) implies that: 

(2p+l)T 2q+1 (-l\iJj 2p+1 r{A) - (2? + l)/ap+i<-l|^g+i^*(-A) = (61) 
Indeed, applying Q to this vector gives 

(-l|^I 2p -l^* 29-1^28+1 Qo = / 2 p+l(-l 1*0* 29-1^28+1 - hq+l (-11*0* 2p-l V>2s+1 = 

Now notice that 

(2,s + l)(-l|V/ 2p _ 1 V2,+i - (2p+l)(-l|C2 S -i^2 P+ i + (-1|VI 2s -iC 2p -iCo 
Hence having in mind the equation (\I/-5| Co = one gets 

(2p+ l)7 2 g+l (-11^+1^*28-1 - (2g+ l)/2p+l(-l|^29+lV'*2 s -l = 

Since it is true for every s we can write it for the generating function as in the eq.(|6l|). 
The second step consists in computing the following average 

A-^-ll^p+iV*^) g{B) «VO)i • • • W))n\ -l-2n) 
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where 



Noticing that 

dDD' 2n+ ^^vm-j==^{D)r(A) 

and calculating the matrix element in a usual way we get the answer: 

A-^-ll^p+i^A) g(B) • • • - 1 - 2n> = A" 1 ^ det(M(A)) 

where M(A) is (n + 1) x (n + 1) matrix with the following matrix elements 

r £)2(i-i) 

M(A) i7 -= / dD. i,j = l,---,n 

v ' J y 

M(A), n +i = — , i = 1, • • • , n 

M(A) n+1 , n+1 - -j== 

Qp04 2 )= y dD w —^y/pm = [VrW)A 2p 

\D\>\A\ 

where [■■•]+ means taking the polynomial part in the expansion around infinity. It is quite obvious that the 
normalized differential dui2 P +i{A) is given by 

du>2 P +i{A) = (2p+ 1)A" 1 det{M(A))dA 

which finishes the proof of the proposition. Returning to the beginning of the proof one find that the 
expression 

d A 

(2 P +i)(-i\^ 2p+1 r(A)— 

can be considered as "symbol" of the normalized differential dui2 P +i(A) . ■ 

The equation {M/ 3 1 Qo = for more complicated states (^-3! than those considered in the Proposition 

16 implies other linear partial differential equations for Bj, so, we get the whole Whitham hierarchy. However 
the equations (|5^) are the only ones with derivatives with respect to only two times. We shall not go further 
into the study of the Whitham hierarchy, because it is not our goal. What we really wanted to do was to 
show the remarkable parallel between the Whitham method and the quantum form factor formulae. We 
hope that this goal is achieved. 

Acknowledgement We would like to thank Tetsuji Miwa for his interest in this work and careful reading 
of the manuscript. 

6 Appendix A 

In this appendix we explain why the conditions 

L%\A 1 ,---,A n \B 1 ,---,B 2n ) 

B 2 »=-S 1 , A n =±B 2n -i 

= —e ± L^ 1 \Ai, ■ • • , A n -x\B 2 , • • • , flan-i) (62) 
(e = + or — respectively for the operators $2fc and their descendents or ^fc+i and their descendents) and 

(n 2n \ 
]J]]^{A t ,B 3 )]l{A 2 -Af)L%\A 1 ,---,A n \B lr --,B 2n )a k n \ = 0, k > n + 1 (63) 
i=l j=l i<j J 
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are sufficient for the locality of the operator whose form factors are given by 

fo(01tfh, ■ ■ • > 02n) + ■■■+ = 

n 2n 

i<j i=lj =n +l yrJ rl ' 
*-Jo{Pl,fa,--,02n) + ■■■ + 

where 

fo(Pl, 02, ■ ■ ■ , 02n) + ■■■+ = 

IP r n 2n n 

= 72-yr / ^ ■ ■ ■ / ^ n n ^(Ai> B i)u.(^ - A f) ■■■,A n \B 1 ,.-., B 2n) n op 

\ > J J i—Xj—x i<j i=l 

The calculations which we are going to make are well known even for more general case ||. However, we 
want to repeat them for our particular situation for the completeness of the exposition. 

In the case of diagonal scattering the only non-trivial requirement for the form factors is the following: 

27ri resp 2n= p 1+ T T ifo{0x, 02, • • • , 02n-l, • ■ - fen) +■■■+ = 



' 2n-l 



3=2 



where the S-matrix in our case is 



Using the identity 



q(R fl s tt sinhl(/3 2 -/? 1 + ^) _ m-bm 

s{P2 ~ Pl) l\ ^\(0 2 -0x-^) - muB 2 ) 



exp(--(i/ -1)(A + &))- 



2 V (sinhi/(^ -/3i)) 2 " ^{B^Bj) 

one find that the relation ( |64| ) is equivalent to 

fo{0X,02, ■ ■ ■ , 02n-X,02n) + ■■■ + 



/ 2n-l 2n-l \ 

= —fo(02,---,02n-i)--+-+ II H-BuB^-e I] V^Bi,^) 

1 \J=2 j=2 / 

if the constant c is taken as c = 2i/(£(— 7ri)) _1 . Explicitly the LHS of this equation is 



n 2n 1 

v ' ^ ^ i=l * 1 1=2 



n(A? - Af) L<g\Ax, ■ ■ -,A n \Bx, ■ ■ -,B 2n _x, -B x ) JJ Q r< (65) 

£<j i= 1 



where we have used the identity 

ii(A,B)^(A,-B)= 



A 2 -B 2 



Let us consider the integral over A n . If the contour is such that \A n \ > \B\\ we can replace in this integral 
a~ n by b^ n+1 (a n - bx)' 1 . Indeed 

,-n+l n-1 
V ^ 3=1 j>n+l 
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n-1 

the sum ^ can be omitted due to anti-symmetry with respect to a,j, j = 1, ■ • • ,n — 1, the sum ^ can 

j=l j>n+l 

be omitted because the integrand does not have residue at A n = oo due to (|63|). Thus the integral over A n 
becomes 

1 f 1 2n_1 

K n+1 y- / dAn -p— p ^(A„,S,)[](4 2 -^) L^>(A 1 ,--- > A n \B 1 ,---,B 2n _ 1 ,-B 1 ) = 

1 /2n-l 

= 2 6 r n+1 Sr 1 II(^ 2 - i? i 2 ) II ^(-fli,5j)4 n) (4-,Vi-Bi|Bi I - 1 52 n -i 1 -5i)- 

i<n V j=2 

2n-l ^ 

- J[ xP(B 1 ,B j )L ( g ) (A u ---,A n - 1> Bi\B 1) --- ) B 2n _ 1 ,-B 1 ) 

3=2 

Let us substitute this expression into (|65|): 

n-1 2n-l 
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k 72^t /w • • / n n n^ 2 - ^r- +i n ^ - *o 

1 ^ ' J J i=l j=2 i<j t=l 

!n-l 

JJ ^(-Bi, B^oW • ' ' . A n-i, -Bi\B x , • ■ • , B 2n _!, 

3=2 

2n-l ^ 

I [ 1 P(Bi,Bj)L^(Ai, • • ■ , A„_i, Bi|Bi, • ■ ■ , B 2n -i, — -Bi) 

1 1 „ „ n— 1 2n— 1 it— a 

(2^r y ^ • • • j Mn-i n n n^ 2 - ^ n * 



X 



i=l j=2 i<j i=l 



2n-l 

-(„) 



J=2 

2n-l \ 

- J| V(Si,B,)i^ ) (A 1 ,---,A„_ 1 ,B 1 |B 1 ,---,B 2n _ 1 ,-B 1 ) (66) 

.=2 / 

where we have replaced b^ 1 a~ t (ai ~ bi) by — a~ l for the following reason: for a\ this expression is b^ 1 — a7_, 
the integrand with b^ 1 is regular at zero, for a% with i > 1 we use anti-symmetry. The final formula (pfl) 
shows that the equation (162) is sufficient for locality. 



7 Appendix B 

The reflectionless case is a rather degenerate one, so, the deformed Riemann bilinear identity |i] does not 
exist in complete form. However for our needs we use only the consequence of the deformed Riemann 
bilinear identity which allows a simple proof in the reflectionless case. We give this proof here for the sake 
of completeness. We have used the following fact: 

. . 2 2n 

dA 1 dA 2 Y[Y[ij(A l ,B J )C(A 1 ,A 2 )a h 1 a l 2 = Vfc,Z (67) 
J J i=ij=i 

where we had the expression for C(A\ 1 A 2 ): 

M) = A^A 2 {iTT^ (P(Al)P(A2) - n-A^m-A^) + (P(-A!)P(A 2 ) - P(A 1 )P(-A 2 ))| (68) 
Let us introduce the functions 

2n 2n 

F(A) = [] ip(A, Bj)P(A), G(A) = JJ ^A, B,)P(-A) 

3 = 1 3=1 
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4>{M,B)= (#-4-1 i>(AB), (Mi) 



Recall that the function VK-^i -B) satisfies the difference equation 

' B — A 
B + qA 

which implies that 

F(Aq) = G(A) 

The integral ( |67|) can be rewritten as follows: 

jit J ^{^^(HAi)F(A 2 )-F( q A 1 )F( q A 2 )) + (F( q A 1 )F(A 2 )-F(A 1 )F( q A 2 ))^ k 1 a l 2 

Changing variables Aj — > gA; where needed one easily find that this integral equals zero. Recall that a, = A? y 
and g 2ly = 1, so a\a l 2 do not change under these changes of variables. 

Similar proof for the case of generic coupling constant is more complicated because one crosses singularities 
when changing variables and moving contours. 

Now let us prove the Proposition 1. We want to find equivalent expressions for C{A\, A 2 ). Let us rewrite 
the expression (|6q) in the integral form: 



C(A 1 ,A 2 ) = C'(A 1 ,A 2 )-C'(A 2 ,A 1 ) 

where 

C'(A A \ = — ( p (M)P{A 2 ) PjA^Pj-A,) P(-A 1 )P(A 2 ) Pj-A^Pj-A,) 
1 u 2> 2A 2 \ A 2 + A\ A 2 -Ax A 2 -Ax A 2 +A 1 

= M f dD 2 f dDl WW 



{D x + D 2 ){Dl-Al){Dl-AD 

\D 2 \>\D 1 \ \Dt\y\At\.\A2\ 



Let us modify this expression by adding "exact forms" in variables A±. It is convenient to use the formula 
QiAjPiAt) - q Q{ q A 1 )P{-A 1 ) = 

= J dD, (Q(J>i)(£>x + A,) - qQ(-qD 1 )(D 1 - A,)) 



(Df-AfY 

\Di\>\A t \ 

Suppose that the polynomial Q(A) solves the equation 

Q(D l)+ qQ(-qD 1 ) = J dD 2 (Di+ ^_ Al) (70) 

|23 a |>|£)i| 

(obviously the RHS of this equation is a polynomial) then we can rewrite the expression for C'(Ai,A 2 ) in 
the following equivalent form 

C'(A U A 2 )= J dD 1 D 1 ^} Dl) A2) (Q(Di) - qQi-qDj)) 

\Di\>\Ax\ 1 

Now we have to solve the equation (|70|). It is simple: 

\D 2 \>\D X \ 

where r](x) satisfy rj(x) — rj(—qx) — j^- Namely, 

w( x ) = ^2 i z a 2k+i x2k ~ ^2 t _ a 2k x2k 

k=0 H fc=l y 



3(i 



Hence 



C'(A 1 ,A 2 ) = J dD 2 J dD 1 P{D 1 )P(D 2 )r e 



|£> 2 |>|£)i| 

where T e (x) = rj(x) + r](—qx), ie. 



\Di\>\A ± \ 



D 2 J (Di-A{)(D'i-Ai) 



1 1 __2fc-l 1 + y 



>( x } 1 + q 2k-i x2k 1 i 



,2/>- 



k=l 



k=l 



The expression for C e {A\, A 2 ) given in Proposition 1 follows from these formulae. The expression for 
C (Ai,A 2 ) can be obtained in a similar way eliminating the even degrees of A\ from C'(A\, A 2 ). 
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